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ABSTRACT 

We evaluate time series observations of chromospheric lines (Mg n, Mg i, and C ii) for the K giant 
a Tau obtained using the WE LWP camera at high dispersion. These observations cover a time span of 
about 2 weeks in 1994 February-March and were designed to resolve variations occurring within hours, 
days, and weeks. We consider the observational results in relation to theoretical acoustic heating models, 
motivated by the fact that a Tau may exhibit a basal (i.e., minimum) level of chromospheric activity. The 
data reveal flux variations between the extremes of 8% in Mg n h + k and 15% in each emission com- 
ponent. These variations occur on timescales as short as 8 hr but not on timescales longer than ~ 3 
days. For the h and k components, flux variations occurring on a timescale as short as 1.5 hr are also 
found. These changes are often not correlated (and are sometimes even anticorrelated), leading to 
remarkable differences in the h/k ratios. We argue that these results are consistent with the presence of 
strong acoustic shocks, which can lead to variable Mg n line emission when only a small number of 
strong shocks are propagating through the atmosphere. We deduce the electron density in the C n] 
A2325 line formation region to be log N e ~ 9.0, in agreement with previous studies. Our data provide 
evidence that the Mg ii basal flux limit for K giants might be a factor of 4 higher than suggested by 
Rutten et al. 

Subject headings: stars: activity — stars: chromospheres — stars: individual (a Tauri) — 
stars: late-type — ultraviolet: stars 


1. INTRODUCTION 

The outer atmospheres of K giant stars have been studied 
extensively with satellite-based instruments over the past 
two decades. Many of these stars are considered to be 
“inactive” as they have relatively low chromospheric emis- 
sion losses and do not possess detectable coronae (e.g., 
Ayres, Fleming, & Schmitt 1991; Ayres et al. 1995). For 
a Tauri (Aldebaran = HD 29139 = HR 1457; K5 III), im- 
portant information about the thermodynamic structure of 
the outer atmosphere has been obtained from WE (Kelch 
et al. 1978; Judge 1986) and Hubble Space Telescope ( HST ) 
Goddard High Resolution Spectrograph (GHRS) data 
(Carpenter et al. 1991). It has been found that the outer 
atmospheres of a Tau and other inactive K giants are domi- 
nated by chromospheres with temperatures between 5000 
and 15,000 K, perhaps accompanied by tiny layers at some- 
what higher temperatures. Stellar winds with moderately 
large mass-loss rates, of order 10” 9 M 0 yr” 1 (or somewhat 
below), are also present. The final flow speeds of the winds 
are relatively low, probably not exceeding 50 km s “ x . Judge 
& Stencel (1991) have evaluated these findings and have 
discussed links to the evolutionary status of the stars by 
interpreting the thermodynamic properties of the atmo- 
spheres and winds by means of statistical concepts. None- 
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theless, the underlying physical processes responsible for 
these properties remain only barely understood, partly 
owing to the lack of theoretical models and suitable 
chromospheric line diagnostics (see, e.g., Jordan & Linsky 
1987). Ideally, chromospheric line diagnostics can give 
observational insight into the height-dependent atmo- 
spheric energy and momentum input, which are key to the 
chromospheric phenomenon. 

Why should we study the outer atmospheric dynamics of 
inactive K giant stars at all? One basic reason is that these 
stars represent the far distant future of solar-type activity. 
When solar-type stars evolve up the giant branch, their 
angular momentum subsides, in part because of the increase 
in the moment of inertia and magnetic braking resulting 
from the onset of massive stellar winds (Schrijver 1993; 
Charbonneau, Schrijver, & MacGregor 1996). Most 
intriguingly, the stellar rotation periods of K giants are 
found to be extremely long, e.g., 223 days for a Boo (K1.5 
III) and 643 days for a Tau (Hatzes & Cochran 1993). 
During the course of post-main-sequence evolution, the 
stellar magnetic fields are expected to weaken dramatically 
and are also likely to change their topology, which is found 
to be relevant to the outer atmospheric dynamics (Rosner et 
al. 1991, 1995). It has been argued that the dominant 
chromospheric heating mechanism in inactive late-type 
giants might be essentially acoustic rather than magnetic. 
This result has found support in studies by Schrijver (1987), 
Pasquifii, Pallavicini, & Pakull (1988), Dupree, Hartmann, 
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& Smith (1990), Rutten et al. (1991), and Strassmeier et al. 
(1994), which show that most single noncoronal giants and 
supergiants have chromospheric fluxes, in agreement with 
or close to the empirically deduced basal flux limits, a 
suspected indicator of pure acoustic heating. Judge & 
Stencel (1991) extrapolated the Mg ii basal flux limits to low 
effective temperatures in order to evaluate the activity level 
of noncoronal K giant stars. These limits were successfully 
reproduced by detailed theoretical ab initio models based 
on acoustic energy dissipation (Cuntz, Rammacher, & 
Ulmschneider 1994; Buchholz & Ulmschneider 1994; 
Buchholz 1995). The agreements were found to be within a 
factor of 2 for stars of widely different spectral types and 
luminosity classes. On the other hand, evidence also exists 
indicating that even inactive noncoronal giants and super- 
giants might have some magnetic field coverage remaining. 
These magnetic fields may be important for explaining the 
observed level of atmospheric turbulence (Judge & Cuntz 
1993) and the stellar mass-loss behavior (Hartmann & 
MacGregor 1980; Cuntz 1990b; Sutmann & Cuntz 1995). 

In order to interpret the chromospheric line emission in 
terms of physical parameters (e.g., energy and momentum 
balance of the chromosphere and wind), studies of line for- 
mation in low-density plasma are ultimately required. 
These studies yield empirical results concerning the forma- 
tion of individual chromospheric lines. For a Tau, results 
were given by Stencel et al. (1981), Brown & Carpenter 
(1984), Carpenter, Brown, & Stencel (1985), Judge (1986), 
and Byrne et al. (1988). The empirical technique gives infor- 
mation about thermodynamic quantities present in the line 
formation region. Nonetheless, this technique is only of 
limited use as the lines in question are usually formed over a 
large range of mass column density (Judge 1990) and are 
furthermore influenced by time-dependent effects. There- 
fore, it is often unclear what the global measures of physical 
properties represent. A recent example demonstrating that 
highly resolved C n] line profiles can provide information 
about chromospheric substructure velocities has been given 
by Judge (1994), who reanalyzed the GHRS data for a Tau. 
Judge found significant evidence for multiple density com- 
ponents in this star moving inward and outward with 
respect to the stellar surface. 

Other studies of stellar outer atmospheres involve the 
calculation of semiempirical models. These models give 
information about the run of the velocity and the various 
thermodynamic quantities in the atmosphere, which are 
assumed to be smooth functions of the atmospheric height. 
These functions are varied in order to fit observed line 
strengths and profiles. In the case of a Tau, Kelch et al. 
(1978) found evidence that the atmospheric temperature 
rises from 2700 to 12,000 K between 3.0 x 10 1 and 
9.5 x 10“ 5 g cm -2 mass column density. However, the rele- 
vance of these models is seriously in question because the 
lines studied (e.g., Ca n H + K) are often formed in a highly 
inhomogeneous environment that contains an irregular dis- 
tribution of high-temperature regions. In this case, the 
chromospheric line formation cannot be assumed to be 
associated with a time-independent "‘mean"’ atmosphere. It 
is indeed possible that the time-averaged temperature 
decreases with increasing atmospheric height (or decreasing 
mass column density), despite the fact that the atmosphere 
contains a significant number of high-temperature spikes 
(see, e.g., Carlsson & Stein 1995), which are prominently 
responsible for reproducing the observed Mg ii and Ca n 


emission. Detailed observations of a Tau and other giants 
show that the atmospheres also contain plasma at both low 
and high temperature, as revealed by observations of the 
fundamental rotation-vibration band of CO at 4.6 //m 
(Wiedemann et al. 1994) and variable absorption and emis- 
sion in He I /10830 (O'Brien & Lambert 1986). Unified 
models for these different types of inhomogeneities are, 
however, clearly beyond the scope of current cool star 
atmospheric models. 

The work presented in this paper represents an attempt 
to constrain the degree of chromospheric variability and 
inhomogeneity in a Tau, as manifested in the variability 
patterns of various chromospheric emission lines. We 
compare the observed line variability with previous theo- 
retical studies and evaluate the overall chromospheric activ- 
ity level. In § 2, we give details of the observations. Section 3 
discusses the data reduction and analysis. In § 4, we give 
results regarding our Mg ii, Mg i, and C ii observations. We 
also comment on existing theoretical models. Section 5 
offers our conclusions. 

2. OBSERVATIONS 

On 12 days, between 1994 February 25 and March 11, 
a Tau was observed with the IUE LWP camera in high 
dispersion. We accumulated a total of 50 images over 6 US 

1 and 8.5 US 2 shifts. The images are divided into three 
categories according to exposure time. We obtained 23 
short exposures of 2.5 -3 minutes (depending on the suspect- 
ed radiation level) to get information about variations in the 
Mg ii line cores and to measure the Mg ii emission fluxes. 
We also obtained 24 images of moderate exposure length 
(10-14 minutes) to search for possible changes in the Mg ii h 
and k line wings. The cores of the Mg n lines in these images 
are overexposed, as expected. Finally, three images with 90 
minute exposure times were taken in order to properly 
expose the C n] >12325 multiplet and the Mg i A2852 line, 
which are a factor of 20 weaker than the Mg ii lines. The 
C ii] multiplet is particularly useful for obtaining informa- 
tion about the chromospheric electron density. 

3. DATA REDUCTION AND ANALYSIS 

The spectra were reduced interactively at the University 
of Colorado, Boulder, from the standard IUE Spectral 
Image Processing System (SIPS) files provided to the Guest 
Observer. All images were processed with the new LWP 
intensity transfer function (LWP ITF2), described by Cassa- 
tella, Lloyd, & Gonzales Riestra (1988). The absolute flux 
calibration was determined using the LWP ITF2 inverse 
sensitivity function S l 1 and the “ new ” high-dispersion cali- 
bration function C(A) (Cassatella et al. 1990). The fluxes 
were corrected both for short-term variations in sensitivity 
caused by changes in camera temperature using the IUE 
SIPS routine TEMPCOR, and for the long-term sensitivity 
degradation of the LWP camera by using LWP images 
obtained after epoch 1984.5 (Garhart 1992). 

We searched for possible variability in the fluxes of the 
Mg ii h and k resonance lines, using only the short-exposure 
LWP spectra in which the Mg ii emission features were 
unsaturated. These observations were spread unevenly over 

2 weeks in an attempt to detect variations on timescales 
ranging from hours to weeks. Considerable care must be 
taken in analyzing the Mg n features if small-scale fluctua- 
tions are to be detected. Following Ayres et al. (1995), we 
minimized the influence of subjective errors on our flux 
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measurements by registering the Mg n emission profiles to a 
common wavelength scale, then integrating the flux over a 
fixed bandpass. The profiles were brought to the common 
wavelength scale through an iterative process. First, the 
initial line profiles were co-added to form a mean profile. 
The individual profiles then were cross-correlated against 
the mean, in order to establish the offset of each profile from 
the reference. In order to avoid interpolation between wave- 
length bins, and the spectral smoothing that results, the 
profiles were shifted by integral numbers of bins. The shifted 
profiles were used to compute a new mean profile, and the 
entire process was repeated until the residual shifts con- 
verged to zero. 

The integration bandpass was chosen by a visual exami- 
nation of the composite profile, constructed by co-adding 
the registered profiles. The bandpass was then applied uni- 
formly to the individual registered profiles to derive the flux 
measurements summarized in Table 1. We derived the 1 a 
error estimates given in Table 1 by assuming that the 
photometric noise appropriate to typical LW-HI exposures 
obeys a simple two-component distribution: a « a 0 
+ <Ti(IJI mdX ) iI2 . The first term refers to the often dominant 
component of the noise that is independent of the local 
intensity and that arises from the background of the 
Vidicon image contributed partly by the camera pedestal 
and partly by radiation “ fogging.'’ The second term refers 
to the component of the photometric error that depends on 
the intensity, analogous to the behavior of a photon- 
counting sensor with negligible background. We estimated 
the two coefficients empirically in each exposure by examin- 
ing the rms deviation of the profile fluxes from those of a 
reference line shape scaled to the same integrated flux. This 
approach assumes that the shape of the line does not 
change. In the presence of significant morphological profile 
variations, the method overestimates the photometric noise. 
In the empirical “calibration” of the noise model, we com- 
pensated for any prior smoothing of the oversampled pro- 


files and for the intrinsic spatial averaging properties of the 
detectors (see discussion by Ayres 1993). 

4. RESULTS AND DISCUSSION 
4.1. Mg II 

The Mg ii h and k lines are particularly useful because 
they are easily accessible diagnostics that have sufficient 
opacity over a wide temperature range (4000 K < T e < 
8000 K) to reveal detailed information on conditions, from 
the low chromosphere through to the outermost detectable 
parts of the wind. The Mg ii h, /c, and h + k fluxes, as well as 
the h/k ratios of a Tau from the 23 short-exposure images, 
are presented in Figure 1. Figure 2 shows, in more detail, 
the series of observations made between JD ( + 2,440,000) 
9415.5 and 9417.5, which has particularly narrow time 
spacing. 

The data give unique insight into the variability of the 
Mg ii chromospheric emission of a Tau. The mean fluxes 
measured for the Mg H lines are F h = 7.37 x 10“ 11 ergs 
cm" 2 s' 1 and F k = 11.38 x 10 -11 ergs cm -2 s -1 . Between 
the extremes, the flux varies by 15% in each emission com- 
ponent and by 8% in Mg ii h + k. The variations are found 
to occur on timescales as short as 8 hr, but not on signifi- 
cantly longer timescales such as several days or weeks. Most 
surprisingly, it is found that the changes in Mg n h and k are 
often not correlated and, at times, appear to be anti- 
correlated. A chromospheric brightening caused by the 
passage of a shock should appear first in the Mg ii h line, 
which, being a factor of 2 less optically thick than the k line, 
is sensitive to conditions deeper in the chromosphere. Flux 
variations observed in images obtained on JD 2,449,417 are 
at least consistent with such a scenario. On this day, an 
increase of 6% or more in the total Mg n flux has occurred, 
enduring at least 8.5-12 hr. This increase is clearly visible in 
both the Mg ii h and k lines, but it is not possible to 

uniquely identify a temporal offset between the changes in 
1 


Mg u A2800 Emission Line Fluxes of a Tauri 


Image 

(LWP) 

DOY 

(1994) 

Exposure 

Start 

(UT) 

Exposure 

Time 

(minutes) 

JD 2,440,000 + 

Flux 

(10 -11 ergs cm -2 s -1 ) 

Flux Ratio 
h/k 

h 

k 

27492 

55 

23:51 

2.5 

9408.494 

7.54 (0.35) 

11.69 (0.41) 

0.645 (0.038) 

27495 

56 

2:26 

2.5 

9408.602 

7.30 (0.11) 

11.34 (0.24) 

0.643 (0.017) 

27505 

56 

19:51 

3.0 

9409.327 

7.77 (0.10) 

11.39 (0.28) 

0.682 (0.019) 

27510 

57 

0:28 

2.5 

9409.520 

7.34 (0.15) 

11.43 (0.34) 

0.642 (0.023) 

27520 

57 

17:53 

3.0 

9410.245 

7.47 (0.29) 

10.88 (0.36) 

0.687 (0.035) 

27526 

57 

23:30 

3.0 

9410.479 

7.63 (0.36) 

11.17 (0.35) 

0.683 (0.039) 

27557 

62 

23:24 

3.0 

9415.475 

7.24 (0.11) 

11.44 (0.18) 

0.633 (0.014) 

27565 

63 

14:48 

3.0 

9416.116 

7.04 (0.17) 

11.29 (0.33) 

0.623 (0.024) 

27577 

64 

1:38 

3.0 

9416.568 

7.53 (0.21) 

11.10 (0.25) 

0.678 (0.024) 

27580 

64 

11:58 

3.0 

9416.999 

6.91 (0.10) 

11.45 (0.19) 

0.604 (0.013) 

27583 

64 

15:49 

3.0 

9417.159 

7.53 (0.09) 

11.18 (0.15) 

0.674 (0.012) 

27585 

64 

17:17 

3.0 

9417.220 

7.35 (0.11) 

11.36 (0.13) 

0.647 (0.012) 

27587 

64 

18:38 

3.0 

9417.276 

7.25 (0.20) 

11.77 (0.10) 

0.616 (0.018) 

27589 

64 

20:03 

3.0 

9417.335 

7.44 (0.12) 

11.66 (0.14) 

0.638 (0.013) 

27591 

64 

21:29 

3.0 

9417.396 

7.95 (0.09) 

11.39 (0.10) 

0.698 (0.010) 

27593 

64 

22:53 

3.0 

9417.453 

7.46 (0.08) 

11.66(0.17) 

0.640 (0.011) 

27595 

65 

0:16 

3.0 

9417.511 

7.67 (0.08) 

11.87 (0.14) 

0.646 (0.010) 

27615 

67 

11:58 

3.0 

9419.999 

7.27 (0.09) 

11.41 (0.13) 

0.637 (0.011) 

27622 

67 

19:57 

3.0 

9420.331 

6.84 (0.11) 

11.16 (0.16) 

0.613 (0.013) 

27632 

68 

20:11 

3.0 

9421.341 

7.18 (0.15) 

11.29 (0.12) 

0.636 (0.015) 

27638 

69 

1:47 

3.0 

9421.574 

7.16(0.13) 

11.27 (0.17) 

0.636 (0.015) 

27645 

69 

20:28 

3.0 

9422.353 

7.37 (0.19) 

11.49(0.16) 

0.641 (0.019) 

27650 

70 

1:04 

3.0 

9422.545 

7.33 (0.12) 

11.15 (0.23) 

0.658 (0.017) 


Note. — The numbers in parenthesis indicate 1 a error estimates. 



9410 9415 9420 

JD(+2440000) 


Fig. 1. — Integrated Mg II fluxes obtained during 2 weeks of IUE observations in 1994 February March. Panel shows the Mg ii h, k, h + k fluxes, and the 
h/k flux ratios. Error bars represent 1 a standard deviations. Note the extremely nonequidistant spacing of the data. The subset of data displayed in Fig. 2 is 
indicated by the window bordered by dashed lines. 



9415 9416 9417 9418 

JD( + 2440000) 


Fig. 2. — Same as Fig. 1, but now for a subset of the data obtained during a 2 day time span 
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these components. In fact, they show uncorrelated flux 
variations with a timescale of 1.5 hr or less. Also, the dura- 
tion of this Mg ii flux increase cannot be precisely deter- 
mined. Because of the data gaps, it can only be estimated 
that this increased emission might persist over a time span 
between 0.5 and 3 days. Figure 3 shows that the flux 
increase is confined to the Mg ii line cores, which form in 
the outermost parts of the line formation region. 

We also found another instance of Mg ii h + k flux 
increase occurring between JD 2,449,420 and JD 2,449,422. 
During that time, a Mg n flux increase of about 5% occurs. 
Because of the lack of data points, however, this increase is 
less significant than that found on JD 2,449,417. In fact, it is 
possible that the increase persisted only 24 hr instead of 48 
hr. This increase is preceded by a rapid decrease in the Mg n 
h + k flux of 4% occurring over a time span of 8 hr. This 
change is particularly well pronounced in the Mg ii h line. 

We searched for variations in the Mg ii h and k line 
wings, using the 24 moderate-length exposures obtained 
concurrently with the shorter exposures previously dis- 
cussed. We visually examined the overplotted Mg ii profiles 
from all 24 images but were unable to find any significant 
changes in the Mg ii line wings. We attribute this result to 
radiative transfer effects. Since photons in the line core have 
first to migrate to the line wings in order to escape, they are 
subjected to multiple scattering. The scatterings will 
average out local variations in the line flux due to the pres- 
ence of different thermodynamic conditions deep within the 
chromosphere. 

A further issue is the overall Mg n activity level. 
Assuming an angular diameter for a Tau of 20.2 mas (Di 
Benedetto & Rabbia 1987), our minimum Mg ii /i + /c flux 


corresponds to a minimum surface flux of 7.7 x 10 4 ergs 
cm 2 s' 1 . This value for the surface flux agrees to within 
10% of that found by Judge (1986). If a Tau is indeed a 
“chromospheric basal flux star” (see § 1), these observations 
seem to indicate that the basal Mg ii emission flux for K 
giants is about a factor of 4 higher than suggested by Rutten 
et al. (1991). 

4.2. Mg i and C ii 

The Mg i A2852 and C ii] A2325 profiles obtained in the 
three 90 minute exposures are shown in Figures 4 and 5, 
respectively. The Julian dates ( + 2,440,000) at the beginning 
of the observations were 9416.007, 9417.056, and 9420.057. 
As discussed by Stencel et al. (1981), the intensities of the 
C ii intersystem multiple! near 2325 A depend on the elec- 
tron density N e in the range 10 9 < N e < 10 7 and may be 
used as electron density diagnostics in cool giants and 
supergiants. Ideally, we wish to compute the intensity 

ratios R { ( = F 2325 38 /^ 2328 . lo)> R 2 ( — ^ 2325 . 38 /^ 2326 . 91 )’ 
and R 3 ( = F 2324 . 67 /F 2326 91 ), where F A(A) is the integrated 
flux of the C 11 ] line with laboratory wavelength 2. Each 
intensity ratio leads to a separate estimate of the electron 
density for each image. Unfortunately, the C 11 ] line at 
2326.91 A is blended with the Fe 11 line at 2327.40 A (Judge 
& Jordan 1991), making the line ratios R 2 and R 3 difficult 
to determine. This problem can only be alleviated with the 
use of high-quality data, as given by HST GHRS 
(Carpenter et al. 1991). 

Therefore, we have calculated only the line intensity ratio 
R v Using the atomic data of Lennon et al. (1985), we find 
that log N e = 9.2, 9.0, and 8.7 for our three images LWP 
27564, 27582, and 27617, respectively. The errors in these 



Fig. 3. — Profiles of Mg 11 h and k emission lines. Upper panel displays the time-averaged profiles calculated from the 23 short exposures in the IUE time 
series. Lower panel displays the profiles obtained at the beginning and end of the series of observations made between JD ( + 2,440,000) 9417.0 and 9417.5 (see 
Fig. 2). 
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2850 2852 2854 2856 

Wavelength (A) 


Fig. 4. — Line profiles for Mg i A2852 obtained from three 90 minute 
IV E exposures (LWP 27564, 27582, and 27617). Julian dates ( + 2,440,000) 
at the beginning of the observations were 9416.007, 9417.056, and 
9420.057, respectively. 

numbers are difficult to estimate because of the poor quality 
of the data and the blend by the (potentially time-variable) 
Fe ii line, which prevent us to use R 2 and R 3 for additional 
electron density determinations. Therefore, we do not infer 
time variability from these data. The values are thus consis- 
tent with the results from previous studies of a Tau: Judge 
(1986) determined that log N e = 9.0 ± 0.3; Byrne et al. 
(1988) found log N e = 9.1 ± 0.2; and Carpenter et al. (1991), 
who calculated all three intensity ratios using HST GHRS 
spectra, deduced an electron density of slightly less than 
log N e = 9.0. Preliminary results of another study, given by 
Stencel et al. (1992), also indicate electron densities of about 
log N e = 9.0, commensurate with the other studies. Thus, 
we conclude that the Mg ii variability does not appear to be 



2324 2326 2328 

Wavelengtn (A) 


Fig. 5. — Same as Fig. 4, but now for the C n] /2325 multiple! 


associated with anomalous conditions in the chromosphere 
of a Tau and is probably a common occurrence. 

4.3. Comments on Theoretical Chromosphere Models 

Although initial chromospheric modeling was performed 
on the basis of time-independent theory (e.g., Hartmann & 
MacGregor 1980), it now appears that chromospheric 
heating and mass loss in giant stars are greatly influenced 
by time-dependent phenomena. Cuntz (1987) studied the 
response of the outer atmosphere of a Boo (K1.5 III) to 
acoustic shocks generated in a “realistic” photospheric 
model (Bohn 1984). He found strong time-dependent epi- 
sodes of momentum and energy deposition that control the 
chromospheric heating. Similar models were produced for 
additional stars and have been further improved. Cuntz 
(1990a) calculated theoretical model chromospheres for a 
Tau based on acoustic shock wave dissipation. These 
models do not consider the formation of distinct chromo- 
spheric lines but try to estimate the loss of radiative energy 
by employing a thin-plasma approximation. Such an 
approach is motivated by the work of Judge (1990), who 
showed that lines such as Mg ii k thermalize in the upper 
photospheres of inactive low-gravity stars. These lines are 
thus effectively thin in the entire chromosphere, and opti- 
cally thin approximations (rather than more sophisticated 
techniques such as radiative transfer with partial 
redistribution) may be used in the calculation of the line 
cooling. The calculations of Cuntz (1990a) solve the hydro- 
dynamic equations fully consistently, with the shocks 
treated as discontinuities. 

An important aspect of these chromospheric models is 
the consideration of shock-shock interactions in determin- 
ing the outer atmospheric dynamics. Shock-shock inter- 
actions are caused by height-dependent effects in the 
propagation of acoustic frequency spectra generated by 
stellar convection (Bohn 1984; Musielak et al. 1994; Kumar 
1994). When acoustic wave modes propagate through the 
stellar atmosphere, shocks of different strengths and, there- 
fore, different speeds are produced at different atmospheric 
heights. The range in shock speed is sufficiently broad to 
allow high-speed shocks to overtake, interact with, and 
engulf slower moving ones. Since an overtaking shock 
acquires the strength of the shock it engulfs, its speed 
increases, so that it can overtake more and more shocks in 
front of it. Consequently, the amount of momentum and 
energy deposition into the atmosphere varies drastically 
with height, depending on the atmospheric dynamics and 
the photospheric frequency spectrum. This scenario pro- 
motes the generation of distinct variability timescales for 
chromospheric lines in the case where only a limited 
number of strong shocks exists simultaneously in the atmo- 
sphere. 

Figure 6 shows the expected distribution of chromo- 
spheric heating events in the chromosphere of a Tau at 
different atmospheric heights, as given by an acoustic 
heating model. Since the Mg ii diagnostics form over 
roughly 4-5 orders of magnitude in mass column density, it 
is possible to compare the observed variability timescale 
with the estimates from the theoretical model. In the case of 
the Mg ii h and k line wings, no variations were found. This 
is different for the Mg ii line cores, which are formed around 
10 4 and 10“ 5 g cm 1 mass column density. Figure 6 pre- 
dicts a Mg ii variability timescale of 2 8 x 10 4 s (5 -20 hr) 
for the case in which the Mg n emission is associated with a 
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Fig. 6. — Histograms showing the number of shocks present in the time-dependent ab initio a Tau model chromosphere of Cuntz (1990a) as a function of 
the elapsed time between the shocks at a given mass column density. The input spectrum of waves was introduced at 2.5 x 10“ 2 gem* 2 , with a mean period 
indicated by the arrow, and with a Gaussian width (FWHM) of 10 4 s. Periods are subject to stochastic changes, but systematically increase with height. Mg n 
lines are formed over regions from 10“ 2 (line wings) to 10* 5 g cm - 2 (line cores) mass column density. 


small number of stochastic shocks. This timescale is similar 
to (or moderately larger than) the period in which the 
maximum of the acoustic frequency spectrum occurs, which 
has been introduced into the atmosphere. 

The heating model used to produce Figure 6 assumes a 
photospheric energy spectrum centered at 3.5 x 10 4 s ( = 10 
hr), consistent with traditional models of acoustic energy 
generation, as given by Bohn (1984). Bohn found that the 
maximum acoustic energy flux occurs at the driving period 
P max = €P co = 4nc z €/yg+, where P co is the acoustic cutoff 
period, c z is the sound speed at the top of the convection 
zone, y is the ratio of the specific heats, and g + is the stellar 
surface gravity. The parameter €, which denotes the position 
of the energy maximum in the acoustic frequency spectrum, 
is largely unconstrained due to uncertainties in the calcu- 
lation of the initial acoustic energy spectrum from the 
model of stellar convection (see, e.g., Ulmschneider, 
Theurer, & Musielak 1996). 

Earlier work of Bohn (1984), as well as recent results of 
Ulmschneider et al. (1996), suggests that € ranges between 
0.10 and 0.30 for late-type giants. For the model given in 
Figure 6, we took € = 0.10. In the case that € = 0.30 is 
adopted, the predicted Mg n variability timescale will 
change to 0.6-2.4 x 10 5 s (15-60 hr). This range comfort- 
ably encompasses the Mg ii variability timescale seen in our 
observations. This prediction gives some confidence that 
the timescales of the observed variations are consistent with 


the presence of acoustic shocks propagating through the 
stellar atmosphere. The theoretically deduced distribution 
of the chromospheric heating events is also influenced by 
uncertainties in For the model given in Figure 6, we 
took log g+ = 1.3 (van Paradijs & Meurs 1974; Judge & 
Stencel 1991), which is uncertain by at least a factor of 2. As 
a consequence, possible P max are expected to spread over 
an even larger range than given above. 

A further possible indication of the presence of shocks lies 
in the behavior of the Mg n h/k ratio. It is found that the 
Mg ii h and k fluxes are usually not correlated. It can be 
expected that the different optical depths of the two lines 
allow the change in emission to be detected first in the Mg ii 
h line while it is still too deep in the chromosphere to be 
seen in the Mg n k line. These flux variations appear to be 
confined to the line core and thus occur in the outermost 
Mg ii emitting regions. Changes in the Mg n h/k ratio could 
thus indicate a shift of energy dissipation with atmospheric 
height, which is fully reconcilable with the behavior found 
in theoretical models of acoustic shock wave propagation. 

5. CONCLUSIONS 

The objective of this WE observing program was to 
obtain new data with carefully arranged time spacing in 
order to examine the propagation of energy and momen- 
tum generated in subphotospheric convection zones 
through the atmosphere of a Tau. This program was aimed 
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at gaining further insight into the mechanism(s) responsible 
for the observed atmospheric structures in cool, inactive 
giant stars. The star a Tau has often been considered a 
standard example of a giant with a noncoronal outer atmo- 
sphere, which is expected to be dominated by nonmagnetic 
heating. Global oscillations (if present) are irrelevant to the 
atmospheric gas dynamics in this type of star (Sutmann & 
Cuntz 1995). The dominant factor in the dynamics and 
heating of such chromospheres should therefore be the pro- 
pagation and dissipation of acoustic waves. 

Our principal results can be summarized as follows: 

1. We obtained Mg ii observations with the IUE LWP 
camera at high dispersion in 1994 February-March. The 
observations show flux variations between the extremes of 
8% in Mg ii h + k and 15% in the individual components. 
These variations occur on timescales as short as 8 hr. In one 
instance, we found a steady 6% increase in the Mg ii h + k 
flux occurring over a time span of 8.5 12 hr. Because of 
gaps in the data, the precise duration of this event cannot be 
determined. It can be said, however, that the Mg ii flux 
increase did not persist for more than ~3 days. In another 
instance, we found evidence for a 5% increase occurring 
over a time span of 12-24 hr. 

2. We found that the changes in Mg ii h and k are often 
not correlated, leading to remarkable variations in the h/k 
ratio. This ratio is found to vary between 0.60 and 0.70. The 
uncorrelated flux variations can occur on timescales as 
short as 1.5 hr. It is possible that the different optical depths 
of the two lines allow a change in emission to be detected in 
the h line while it is still too deep in the chromosphere to be 
seen in the k line. The flux variations appear to be confined 
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to the line core. Thus, they occur in the outermost Mg n 
emitting regions and might indicate spatial shifts in the 
atmospheric heating rate. 

3. We deduced the electron density in the C ii] 22325 line 
formation region at the epoch of our observations. We 
found electron densities of log N e ~ 9.0, in agreement with 
results from Judge (1986), Byrne et al. (1988), and Carpenter 
et al. (1991). 

4. We evaluated the level of chromospheric activity in 
a Tau based on the total Mg ii emission flux. We found that 
the Mg ii surface flux agrees with those of other inactive K 
giants (see, e.g., Judge & Stencel 1991). Our data indicate 
that the Mg ii basal flux limit in K giants might be a factor 
of 4 higher than suggested by Rutten et al. (1991). 

5. We argue that these results are consistent with the 
presence of strong acoustic shocks produced by the action 
of convective motions as the variability timescale of the 
Mg ii flux is consistent with those expected from acoustic 
shock wave models. The shocks can lead to variable Mg ii 
line emission when only a limited number of strong shocks 
exist simultaneously in the atmosphere and can account for 
the observed variations in the Mg n h/k ratio. 
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